original/building all combinations of msms.R
In HealthEconomicsHackathon/ceatutorial: What the Package Does (One Line, Title Case)
##########################################################################################
##========================================================================================
## CALCULATING STATE OCCUPANCY PROBABILITIES
##========================================================================================
##########################################################################################
### DUE TO COMPUTATIONAL DIFFICULTIES USING GOMPERTZ AND GENERALISED GAMMA FOR THE FIRST
### TRANSITION, MODELS THAT INCLUDE THIS USE MORE TIME POINTS. ALL OTHER MODELS USE THE
### EQUIVALENT OF MONTHLY CYCLES.
##############################
### RFC TREATMENT ARM
#############################
##############################
### exp3
#############################
ptm <- proc.time() #start counting the time it takes
##### exp2exp3
exp1exp2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "exp", "exp"),
trans=tmat2, M=5000)
(proc.time() - ptm) #stop counting the time it took (seconds)
(proc.time() - ptm)/60 #stop counting the time it took (minutes)
(proc.time() - ptm)/3600 #stop counting the time it took (hours)
wei1exp2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "exp", "exp"),
trans=tmat2, M=5000)
gom1exp2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
dist=cbind("gom", "exp", "exp"),
trans=tmat2, M=5000)
gam1exp2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "exp", "exp"),
trans=tmat2, M=5000)
logn1exp2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "exp", "exp"),
trans=tmat2, M=5000)
logl1exp2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "exp", "exp"),
trans=tmat2, M=5000)
##### wei2exp3
exp1wei2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "wei", "exp"),
trans=tmat2, M=5000)
wei1wei2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "wei", "exp"),
trans=tmat2, M=5000)
gom1wei2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
dist=cbind("gom", "wei", "exp"),
trans=tmat2, M=5000)
gam1wei2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "wei", "exp"),
trans=tmat2, M=5000)
logn1wei2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "wei", "exp"),
trans=tmat2, M=5000)
logl1wei2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "wei", "exp"),
trans=tmat2, M=5000)
##### gom2exp3
exp1gom2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "gom", "exp"),
trans=tmat2, M=5000)
wei1gom2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "gom", "exp"),
trans=tmat2, M=5000)
gom1gom2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
dist=cbind("gom", "gom", "exp"),
trans=tmat2, M=5000)
gam1gom2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gom", "exp"),
trans=tmat2, M=5000)
logn1gom2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "gom", "exp"),
trans=tmat2, M=5000)
logl1gom2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "gom", "exp"),
trans=tmat2, M=5000)
##### gam2exp3
exp1gam2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "gam", "exp"),
trans=tmat2, M=5000)
wei1gam2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "gam", "exp"),
trans=tmat2, M=5000)
gom1gam2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
dist=cbind("gom", "gam", "exp"),
trans=tmat2, M=5000)
gam1gam2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gam", "exp"),
trans=tmat2, M=5000)
logn1gam2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "gam", "exp"),
trans=tmat2, M=5000)
logl1gam2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "gam", "exp"),
trans=tmat2, M=5000)
##### logn2exp3
exp1logn2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "logn", "exp"),
trans=tmat2, M=5000)
wei1logn2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "logn", "exp"),
trans=tmat2, M=5000)
gom1logn2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
dist=cbind("gom", "logn", "exp"),
trans=tmat2, M=5000)
gam1logn2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logn", "exp"),
trans=tmat2, M=5000)
logn1logn2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "logn", "exp"),
trans=tmat2, M=5000)
logl1logn2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "logn", "exp"),
trans=tmat2, M=5000)
##### logl2exp3
exp1logl2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "logl", "exp"),
trans=tmat2, M=5000)
wei1logl2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "logl", "exp"),
trans=tmat2, M=5000)
gom1logl2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
dist=cbind("gom", "logl", "exp"),
trans=tmat2, M=5000)
gam1logl2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logl", "exp"),
trans=tmat2, M=5000)
logn1logl2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "logl", "exp"),
trans=tmat2, M=5000)
logl1logl2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "logl", "exp"),
trans=tmat2, M=5000)
##############################
### wei3
#############################
##### exp2wei3
exp1exp2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "exp", "wei"),
trans=tmat2, M=5000)
wei1exp2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "exp", "wei"),
trans=tmat2, M=5000)
gom1exp2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("gom", "exp", "wei"),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
trans=tmat2, M=5000)
gam1exp2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "exp", "wei"),
trans=tmat2, M=5000)
logn1exp2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "exp", "wei"),
trans=tmat2, M=5000)
logl1exp2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "exp", "wei"),
trans=tmat2, M=5000)
##### wei2wei3
exp1wei2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "wei", "wei"),
trans=tmat2, M=5000)
wei1wei2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "wei", "wei"),
trans=tmat2, M=5000)
gom1wei2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
dist=cbind("gom", "wei", "wei"),
trans=tmat2, M=5000)
gam1wei2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "wei", "wei"),
trans=tmat2, M=5000)
logn1wei2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "wei", "wei"),
trans=tmat2, M=5000)
logl1wei2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "wei", "wei"),
trans=tmat2, M=5000)
##### gom2wei3
exp1gom2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "gom", "wei"),
trans=tmat2, M=5000)
wei1gom2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "gom", "wei"),
trans=tmat2, M=5000)
gom1gom2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
dist=cbind("gom", "gom", "wei"),
trans=tmat2, M=5000)
gam1gom2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gom", "wei"),
trans=tmat2, M=5000)
logn1gom2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "gom", "wei"),
trans=tmat2, M=5000)
logl1gom2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "gom", "wei"),
trans=tmat2, M=5000)
##### gam2wei3
exp1gam2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "gam", "wei"),
trans=tmat2, M=5000)
wei1gam2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "gam", "wei"),
trans=tmat2, M=5000)
gom1gam2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
dist=cbind("gom", "gam", "wei"),
trans=tmat2, M=5000)
gam1gam2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gam", "wei"),
trans=tmat2, M=5000)
logn1gam2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "gam", "wei"),
trans=tmat2, M=5000)
logl1gam2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "gam", "wei"),
trans=tmat2, M=5000)
##### logn2wei3
exp1logn2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "logn", "wei"),
trans=tmat2, M=5000)
wei1logn2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "logn", "wei"),
trans=tmat2, M=5000)
gom1logn2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
dist=cbind("gom", "logn", "wei"),
trans=tmat2, M=5000)
gam1logn2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logn", "wei"),
trans=tmat2, M=5000)
logn1logn2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "logn", "wei"),
trans=tmat2, M=5000)
logl1logn2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "logn", "wei"),
trans=tmat2, M=5000)
##### logl2wei3
exp1logl2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "logl", "wei"),
trans=tmat2, M=5000)
wei1logl2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "logl", "wei"),
trans=tmat2, M=5000)
gom1logl2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
dist=cbind("gom", "logl", "wei"),
trans=tmat2, M=5000)
gam1logl2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logl", "wei"),
trans=tmat2, M=5000)
logn1logl2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "logl", "wei"),
trans=tmat2, M=5000)
logl1logl2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "logl", "wei"),
trans=tmat2, M=5000)
##############################
### gom3
#############################
##### exp2gom3
exp1exp2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "exp", "gom"),
trans=tmat2, M=5000)
wei1exp2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "exp", "gom"),
trans=tmat2, M=5000)
gom1exp2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
dist=cbind("gom", "exp", "gom"),
trans=tmat2, M=5000)
gam1exp2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "exp", "gom"),
trans=tmat2, M=5000)
logn1exp2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "exp", "gom"),
trans=tmat2, M=5000)
logl1exp2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "exp", "gom"),
trans=tmat2, M=5000)
##### wei2gom3
exp1wei2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "wei", "gom"),
trans=tmat2, M=5000)
wei1wei2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "wei", "gom"),
trans=tmat2, M=5000)
gom1wei2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
dist=cbind("gom", "wei", "gom"),
trans=tmat2, M=5000)
gam1wei2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "wei", "gom"),
trans=tmat2, M=5000)
logn1wei2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "wei", "gom"),
trans=tmat2, M=5000)
logl1wei2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "wei", "gom"),
trans=tmat2, M=5000)
##### gom2gom3
exp1gom2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "gom", "gom"),
trans=tmat2, M=5000)
wei1gom2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "gom", "gom"),
trans=tmat2, M=5000)
gom1gom2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
dist=cbind("gom", "gom", "gom"),
trans=tmat2, M=5000)
gam1gom2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gom", "gom"),
trans=tmat2, M=5000)
logn1gom2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "gom", "gom"),
trans=tmat2, M=5000)
logl1gom2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "gom", "gom"),
trans=tmat2, M=5000)
##### gam2gom3
exp1gam2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "gam", "gom"),
trans=tmat2, M=5000)
wei1gam2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "gam", "gom"),
trans=tmat2, M=5000)
gom1gam2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
dist=cbind("gom", "gam", "gom"),
trans=tmat2, M=5000)
gam1gam2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gam", "gom"),
trans=tmat2, M=5000)
logn1gam2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "gam", "gom"),
trans=tmat2, M=5000)
logl1gam2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "gam", "gom"),
trans=tmat2, M=5000)
##### logn2gom3
exp1logn2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "logn", "gom"),
trans=tmat2, M=5000)
wei1logn2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "logn", "gom"),
trans=tmat2, M=5000)
gom1logn2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
dist=cbind("gom", "logn", "gom"),
trans=tmat2, M=5000)
gam1logn2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logn", "gom"),
trans=tmat2, M=5000)
logn1logn2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "logn", "gom"),
trans=tmat2, M=5000)
logl1logn2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "logn", "gom"),
trans=tmat2, M=5000)
##### logl2gom3
exp1logl2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "logl", "gom"),
trans=tmat2, M=5000)
wei1logl2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "logl", "gom"),
trans=tmat2, M=5000)
gom1logl2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
dist=cbind("gom", "logl", "gom"),
trans=tmat2, M=5000)
gam1logl2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logl", "gom"),
trans=tmat2, M=5000)
logn1logl2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "logl", "gom"),
trans=tmat2, M=5000)
logl1logl2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "logl", "gom"),
trans=tmat2, M=5000)
##############################
### gam3
#############################
##### exp2gam3
exp1exp2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "exp", "gam"),
trans=tmat2, M=5000)
wei1exp2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "exp", "gam"),
trans=tmat2, M=5000)
gom1exp2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("gom", "exp", "gam"),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
trans=tmat2, M=5000)
gam1exp2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "exp", "gam"),
trans=tmat2, M=5000)
logn1exp2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "exp", "gam"),
trans=tmat2, M=5000)
logl1exp2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "exp", "gam"),
trans=tmat2, M=5000)
##### wei2gam3
exp1wei2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "wei", "gam"),
trans=tmat2, M=5000)
wei1wei2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "wei", "gam"),
trans=tmat2, M=5000)
gom1wei2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("gom", "wei", "gam"),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
trans=tmat2, M=5000)
gam1wei2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "wei", "gam"),
trans=tmat2, M=5000)
logn1wei2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "wei", "gam"),
trans=tmat2, M=5000)
logl1wei2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "wei", "gam"),
trans=tmat2, M=5000)
##### gom2gam3
exp1gom2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "gom", "gam"),
trans=tmat2, M=5000)
wei1gom2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "gom", "gam"),
trans=tmat2, M=5000)
gom1gom2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("gom", "gom", "gam"),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
trans=tmat2, M=5000)
gam1gom2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gom", "gam"),
trans=tmat2, M=5000)
logn1gom2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "gom", "gam"),
trans=tmat2, M=5000)
logl1gom2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "gom", "gam"),
trans=tmat2, M=5000)
##### gam2gam3
exp1gam2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "gam", "gam"),
trans=tmat2, M=5000)
wei1gam2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "gam", "gam"),
trans=tmat2, M=5000)
gom1gam2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("gom", "gam", "gam"),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
trans=tmat2, M=5000)
gam1gam2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gam", "gam"),
trans=tmat2, M=5000)
logn1gam2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "gam", "gam"),
trans=tmat2, M=5000)
logl1gam2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "gam", "gam"),
trans=tmat2, M=5000)
##### logn2gam3
exp1logn2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "logn", "gam"),
trans=tmat2, M=5000)
wei1logn2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "logn", "gam"),
trans=tmat2, M=5000)
gom1logn2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("gom", "logn", "gam"),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
trans=tmat2, M=5000)
gam1logn2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logn", "gam"),
trans=tmat2, M=5000)
logn1logn2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "logn", "gam"),
trans=tmat2, M=5000)
logl1logn2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "logn", "gam"),
trans=tmat2, M=5000)
##### logl2gam3
exp1logl2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "logl", "gam"),
trans=tmat2, M=5000)
wei1logl2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "logl", "gam"),
trans=tmat2, M=5000)
gom1logl2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("gom", "logl", "gam"),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
trans=tmat2, M=5000)
gam1logl2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logl", "gam"),
trans=tmat2, M=5000)
logn1logl2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "logl", "gam"),
trans=tmat2, M=5000)
logl1logl2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "logl", "gam"),
trans=tmat2, M=5000)
##############################
### logl3
#############################
##### exp2logl3
exp1exp2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "exp", "logl"),
trans=tmat2, M=5000)
wei1exp2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "exp", "logl"),
trans=tmat2, M=5000)
gom1exp2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("gom", "exp", "logl"),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
trans=tmat2, M=5000)
gam1exp2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "exp", "logl"),
trans=tmat2, M=5000)
logn1exp2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "exp", "logl"),
trans=tmat2, M=5000)
logl1exp2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "exp", "logl"),
trans=tmat2, M=5000)
##### wei2logl3
exp1wei2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "wei", "logl"),
trans=tmat2, M=5000)
wei1wei2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "wei", "logl"),
trans=tmat2, M=5000)
gom1wei2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("gom", "wei", "logl"),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
trans=tmat2, M=5000)
gam1wei2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "wei", "logl"),
trans=tmat2, M=5000)
logn1wei2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "wei", "logl"),
trans=tmat2, M=5000)
logl1wei2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "wei", "logl"),
trans=tmat2, M=5000)
##### gom2logl3
exp1gom2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "gom", "logl"),
trans=tmat2, M=5000)
wei1gom2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "gom", "logl"),
trans=tmat2, M=5000)
gom1gom2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("gom", "gom", "logl"),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
trans=tmat2, M=5000)
gam1gom2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gom", "logl"),
trans=tmat2, M=5000)
logn1gom2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "gom", "logl"),
trans=tmat2, M=5000)
logl1gom2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "gom", "logl"),
trans=tmat2, M=5000)
##### gam2logl3
exp1gam2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "gam", "logl"),
trans=tmat2, M=5000)
wei1gam2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "gam", "logl"),
trans=tmat2, M=5000)
gom1gam2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("gom", "gam", "logl"),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
trans=tmat2, M=5000)
gam1gam2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gam", "logl"),
trans=tmat2, M=5000)
logn1gam2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "gam", "logl"),
trans=tmat2, M=5000)
logl1gam2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "gam", "logl"),
trans=tmat2, M=5000)
##### logn2logl3
exp1logn2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "logn", "logl"),
trans=tmat2, M=5000)
wei1logn2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "logn", "logl"),
trans=tmat2, M=5000)
gom1logn2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("gom", "logn", "logl"),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
trans=tmat2, M=5000)
gam1logn2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logn", "logl"),
trans=tmat2, M=5000)
logn1logn2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "logn", "logl"),
trans=tmat2, M=5000)
logl1logn2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "logn", "logl"),
trans=tmat2, M=5000)
##### logl2logl3
exp1logl2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "logl", "logl"),
trans=tmat2, M=5000)
wei1logl2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "logl", "logl"),
trans=tmat2, M=5000)
gom1logl2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("gom", "logl", "logl"),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
trans=tmat2, M=5000)
gam1logl2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logl", "logl"),
trans=tmat2, M=5000)
logn1logl2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "logl", "logl"),
trans=tmat2, M=5000)
logl1logl2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "logl", "logl"),
trans=tmat2, M=5000)
##############################
### logn3
#############################
##### exp2logn3
exp1exp2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "exp", "logn"),
trans=tmat2, M=5000)
wei1exp2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "exp", "logn"),
trans=tmat2, M=5000)
gom1exp2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("gom", "exp", "logn"),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
trans=tmat2, M=5000)
gam1exp2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "exp", "logn"),
trans=tmat2, M=5000)
logn1exp2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "exp", "logn"),
trans=tmat2, M=5000)
logl1exp2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "exp", "logn"),
trans=tmat2, M=5000)
##### wei2logn3
exp1wei2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "wei", "logn"),
trans=tmat2, M=5000)
wei1wei2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "wei", "logn"),
trans=tmat2, M=5000)
gom1wei2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("gom", "wei", "logn"),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
trans=tmat2, M=5000)
gam1wei2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "wei", "logn"),
trans=tmat2, M=5000)
logn1wei2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "wei", "logn"),
trans=tmat2, M=5000)
logl1wei2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "wei", "logn"),
trans=tmat2, M=5000)
##### gom2logn3
exp1gom2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "gom", "logn"),
trans=tmat2, M=5000)
wei1gom2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "gom", "logn"),
trans=tmat2, M=5000)
gom1gom2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("gom", "gom", "logn"),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
trans=tmat2, M=5000)
gam1gom2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gom", "logn"),
trans=tmat2, M=5000)
logn1gom2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "gom", "logn"),
trans=tmat2, M=5000)
logl1gom2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "gom", "logn"),
trans=tmat2, M=5000)
##### gam2logn3
exp1gam2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "gam", "logn"),
trans=tmat2, M=5000)
wei1gam2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "gam", "logn"),
trans=tmat2, M=5000)
gom1gam2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("gom", "gam", "logn"),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
trans=tmat2, M=5000)
gam1gam2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gam", "logn"),
trans=tmat2, M=5000)
logn1gam2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "gam", "logn"),
trans=tmat2, M=5000)
logl1gam2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "gam", "logn"),
trans=tmat2, M=5000)
##### logn2logn3
exp1logn2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "logn", "logn"),
trans=tmat2, M=5000)
wei1logn2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "logn", "logn"),
trans=tmat2, M=5000)
gom1logn2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("gom", "logn", "logn"),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
trans=tmat2, M=5000)
gam1logn2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logn", "logn"),
trans=tmat2, M=5000)
logn1logn2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "logn", "logn"),
trans=tmat2, M=5000)
logl1logn2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "logn", "logn"),
trans=tmat2, M=5000)
##### logl2logn3
exp1logl2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "logl", "logn"),
trans=tmat2, M=5000)
wei1logl2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "logl", "logn"),
trans=tmat2, M=5000)
gom1logl2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("gom", "logl", "logn"),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
trans=tmat2, M=5000)
gam1logl2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logl", "logn"),
trans=tmat2, M=5000)
logn1logl2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "logl", "logn"),
trans=tmat2, M=5000)
logl1logl2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "logl", "logn"),
trans=tmat2, M=5000)
##############################
# ### FC TREATMENT ARM
# #############################
#
# ##############################
# ### exp3
# #############################
#
##### exp2exp3
exp1exp2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "exp", "exp"),
trans=tmat2, M=5000)
wei1exp2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "exp", "exp"),
trans=tmat2, M=5000)
gom1exp2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "exp", "exp"),
trans=tmat2, M=5000)
gam1exp2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)),
dist=cbind("gam", "exp", "exp"),
trans=tmat2, M=5000)
logn1exp2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "exp", "exp"),
trans=tmat2, M=5000)
logl1exp2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "exp", "exp"),
trans=tmat2, M=5000)
##### wei2exp3
exp1wei2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "wei", "exp"),
trans=tmat2, M=5000)
wei1wei2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "wei", "exp"),
trans=tmat2, M=5000)
gom1wei2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "wei", "exp"),
trans=tmat2, M=5000)
gam1wei2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "wei", "exp"),
trans=tmat2, M=5000)
logn1wei2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "wei", "exp"),
trans=tmat2, M=5000)
logl1wei2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "wei", "exp"),
trans=tmat2, M=5000)
##### gom2exp3
exp1gom2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "gom", "exp"),
trans=tmat2, M=5000)
wei1gom2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "gom", "exp"),
trans=tmat2, M=5000)
gom1gom2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "gom", "exp"),
trans=tmat2, M=5000)
gam1gom2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gom", "exp"),
trans=tmat2, M=5000)
logn1gom2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "gom", "exp"),
trans=tmat2, M=5000)
logl1gom2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "gom", "exp"),
trans=tmat2, M=5000)
##### gam2exp3
exp1gam2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "gam", "exp"),
trans=tmat2, M=5000)
wei1gam2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "gam", "exp"),
trans=tmat2, M=5000)
gom1gam2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "gam", "exp"),
trans=tmat2, M=5000)
gam1gam2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gam", "exp"),
trans=tmat2, M=5000)
logn1gam2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "gam", "exp"),
trans=tmat2, M=5000)
logl1gam2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "gam", "exp"),
trans=tmat2, M=5000)
##### logn2exp3
exp1logn2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "logn", "exp"),
trans=tmat2, M=5000)
wei1logn2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "logn", "exp"),
trans=tmat2, M=5000)
gom1logn2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "logn", "exp"),
trans=tmat2, M=5000)
gam1logn2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logn", "exp"),
trans=tmat2, M=5000)
logn1logn2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "logn", "exp"),
trans=tmat2, M=5000)
logl1logn2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "logn", "exp"),
trans=tmat2, M=5000)
##### logl2exp3
exp1logl2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "logl", "exp"),
trans=tmat2, M=5000)
wei1logl2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "logl", "exp"),
trans=tmat2, M=5000)
gom1logl2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "logl", "exp"),
trans=tmat2, M=5000)
gam1logl2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logl", "exp"),
trans=tmat2, M=5000)
logn1logl2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "logl", "exp"),
trans=tmat2, M=5000)
logl1logl2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "logl", "exp"),
trans=tmat2, M=5000)
##############################
### wei3
#############################
##### exp2wei3
exp1exp2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "exp", "wei"),
trans=tmat2, M=5000)
wei1exp2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "exp", "wei"),
trans=tmat2, M=5000)
gom1exp2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "exp", "wei"),
trans=tmat2, M=5000)
gam1exp2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "exp", "wei"),
trans=tmat2, M=5000)
logn1exp2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "exp", "wei"),
trans=tmat2, M=5000)
logl1exp2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "exp", "wei"),
trans=tmat2, M=5000)
##### wei2wei3
exp1wei2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "wei", "wei"),
trans=tmat2, M=5000)
wei1wei2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "wei", "wei"),
trans=tmat2, M=5000)
gom1wei2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "wei", "wei"),
trans=tmat2, M=5000)
gam1wei2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "wei", "wei"),
trans=tmat2, M=5000)
logn1wei2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "wei", "wei"),
trans=tmat2, M=5000)
logl1wei2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "wei", "wei"),
trans=tmat2, M=5000)
##### gom2wei3
exp1gom2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "gom", "wei"),
trans=tmat2, M=5000)
wei1gom2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "gom", "wei"),
trans=tmat2, M=5000)
gom1gom2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "gom", "wei"),
trans=tmat2, M=5000)
gam1gom2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gom", "wei"),
trans=tmat2, M=5000)
logn1gom2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "gom", "wei"),
trans=tmat2, M=5000)
logl1gom2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "gom", "wei"),
trans=tmat2, M=5000)
##### gam2wei3
exp1gam2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "gam", "wei"),
trans=tmat2, M=5000)
wei1gam2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "gam", "wei"),
trans=tmat2, M=5000)
gom1gam2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "gam", "wei"),
trans=tmat2, M=5000)
gam1gam2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gam", "wei"),
trans=tmat2, M=5000)
logn1gam2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "gam", "wei"),
trans=tmat2, M=5000)
logl1gam2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "gam", "wei"),
trans=tmat2, M=5000)
##### logn2wei3
exp1logn2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "logn", "wei"),
trans=tmat2, M=5000)
wei1logn2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "logn", "wei"),
trans=tmat2, M=5000)
gom1logn2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "logn", "wei"),
trans=tmat2, M=5000)
gam1logn2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logn", "wei"),
trans=tmat2, M=5000)
logn1logn2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "logn", "wei"),
trans=tmat2, M=5000)
logl1logn2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "logn", "wei"),
trans=tmat2, M=5000)
##### logl2wei3
exp1logl2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "logl", "wei"),
trans=tmat2, M=5000)
wei1logl2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "logl", "wei"),
trans=tmat2, M=5000)
gom1logl2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "logl", "wei"),
trans=tmat2, M=5000)
gam1logl2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logl", "wei"),
trans=tmat2, M=5000)
logn1logl2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "logl", "wei"),
trans=tmat2, M=5000)
logl1logl2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "logl", "wei"),
trans=tmat2, M=5000)
##############################
### gom3
#############################
##### exp2gom3
exp1exp2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "exp", "gom"),
trans=tmat2, M=5000)
wei1exp2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "exp", "gom"),
trans=tmat2, M=5000)
gom1exp2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "exp", "gom"),
trans=tmat2, M=5000)
gam1exp2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "exp", "gom"),
trans=tmat2, M=5000)
logn1exp2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "exp", "gom"),
trans=tmat2, M=5000)
logl1exp2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "exp", "gom"),
trans=tmat2, M=5000)
##### wei2gom3
exp1wei2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "wei", "gom"),
trans=tmat2, M=5000)
wei1wei2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "wei", "gom"),
trans=tmat2, M=5000)
gom1wei2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "wei", "gom"),
trans=tmat2, M=5000)
gam1wei2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "wei", "gom"),
trans=tmat2, M=5000)
logn1wei2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "wei", "gom"),
trans=tmat2, M=5000)
logl1wei2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "wei", "gom"),
trans=tmat2, M=5000)
##### gom2gom3
exp1gom2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "gom", "gom"),
trans=tmat2, M=5000)
wei1gom2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "gom", "gom"),
trans=tmat2, M=5000)
gom1gom2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "gom", "gom"),
trans=tmat2, M=5000)
gam1gom2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gom", "gom"),
trans=tmat2, M=5000)
logn1gom2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "gom", "gom"),
trans=tmat2, M=5000)
logl1gom2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "gom", "gom"),
trans=tmat2, M=5000)
##### gam2gom3
exp1gam2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "gam", "gom"),
trans=tmat2, M=5000)
wei1gam2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "gam", "gom"),
trans=tmat2, M=5000)
gom1gam2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "gam", "gom"),
trans=tmat2, M=5000)
gam1gam2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gam", "gom"),
trans=tmat2, M=5000)
logn1gam2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "gam", "gom"),
trans=tmat2, M=5000)
logl1gam2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "gam", "gom"),
trans=tmat2, M=5000)
##### logn2gom3
exp1logn2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "logn", "gom"),
trans=tmat2, M=5000)
wei1logn2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "logn", "gom"),
trans=tmat2, M=5000)
gom1logn2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "logn", "gom"),
trans=tmat2, M=5000)
gam1logn2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logn", "gom"),
trans=tmat2, M=5000)
logn1logn2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "logn", "gom"),
trans=tmat2, M=5000)
logl1logn2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "logn", "gom"),
trans=tmat2, M=5000)
##### logl2gom3
exp1logl2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "logl", "gom"),
trans=tmat2, M=5000)
wei1logl2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "logl", "gom"),
trans=tmat2, M=5000)
gom1logl2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "logl", "gom"),
trans=tmat2, M=5000)
gam1logl2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logl", "gom"),
trans=tmat2, M=5000)
logn1logl2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "logl", "gom"),
trans=tmat2, M=5000)
logl1logl2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "logl", "gom"),
trans=tmat2, M=5000)
##############################
### gam3
#############################
##### exp2gam3
exp1exp2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "exp", "gam"),
trans=tmat2, M=5000)
wei1exp2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "exp", "gam"),
trans=tmat2, M=5000)
gom1exp2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "exp", "gam"),
trans=tmat2, M=5000)
gam1exp2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "exp", "gam"),
trans=tmat2, M=5000)
logn1exp2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "exp", "gam"),
trans=tmat2, M=5000)
logl1exp2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "exp", "gam"),
trans=tmat2, M=5000)
##### wei2gam3
exp1wei2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "wei", "gam"),
trans=tmat2, M=5000)
wei1wei2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "wei", "gam"),
trans=tmat2, M=5000)
gom1wei2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "wei", "gam"),
trans=tmat2, M=5000)
gam1wei2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "wei", "gam"),
trans=tmat2, M=5000)
logn1wei2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "wei", "gam"),
trans=tmat2, M=5000)
logl1wei2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "wei", "gam"),
trans=tmat2, M=5000)
##### gom2gam3
exp1gom2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "gom", "gam"),
trans=tmat2, M=5000)
wei1gom2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "gom", "gam"),
trans=tmat2, M=5000)
gom1gom2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "gom", "gam"),
trans=tmat2, M=5000)
gam1gom2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gom", "gam"),
trans=tmat2, M=5000)
logn1gom2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "gom", "gam"),
trans=tmat2, M=5000)
logl1gom2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "gom", "gam"),
trans=tmat2, M=5000)
##### gam2gam3
exp1gam2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "gam", "gam"),
trans=tmat2, M=5000)
wei1gam2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "gam", "gam"),
trans=tmat2, M=5000)
gom1gam2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "gam", "gam"),
trans=tmat2, M=5000)
gam1gam2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gam", "gam"),
trans=tmat2, M=5000)
logn1gam2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "gam", "gam"),
trans=tmat2, M=5000)
logl1gam2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "gam", "gam"),
trans=tmat2, M=5000)
##### logn2gam3
exp1logn2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "logn", "gam"),
trans=tmat2, M=5000)
wei1logn2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "logn", "gam"),
trans=tmat2, M=5000)
gom1logn2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "logn", "gam"),
trans=tmat2, M=5000)
gam1logn2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logn", "gam"),
trans=tmat2, M=5000)
logn1logn2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "logn", "gam"),
trans=tmat2, M=5000)
logl1logn2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "logn", "gam"),
trans=tmat2, M=5000)
##### logl2gam3
exp1logl2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "logl", "gam"),
trans=tmat2, M=5000)
wei1logl2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "logl", "gam"),
trans=tmat2, M=5000)
gom1logl2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "logl", "gam"),
trans=tmat2, M=5000)
gam1logl2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logl", "gam"),
trans=tmat2, M=5000)
logn1logl2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "logl", "gam"),
trans=tmat2, M=5000)
logl1logl2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "logl", "gam"),
trans=tmat2, M=5000)
##############################
### logl3
#############################
##### exp2logl3
exp1exp2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "exp", "logl"),
trans=tmat2, M=5000)
wei1exp2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "exp", "logl"),
trans=tmat2, M=5000)
gom1exp2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "exp", "logl"),
trans=tmat2, M=5000)
gam1exp2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "exp", "logl"),
trans=tmat2, M=5000)
logn1exp2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "exp", "logl"),
trans=tmat2, M=5000)
logl1exp2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "exp", "logl"),
trans=tmat2, M=5000)
##### wei2logl3
exp1wei2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "wei", "logl"),
trans=tmat2, M=5000)
wei1wei2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "wei", "logl"),
trans=tmat2, M=5000)
gom1wei2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "wei", "logl"),
trans=tmat2, M=5000)
gam1wei2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "wei", "logl"),
trans=tmat2, M=5000)
logn1wei2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "wei", "logl"),
trans=tmat2, M=5000)
logl1wei2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "wei", "logl"),
trans=tmat2, M=5000)
##### gom2logl3
exp1gom2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "gom", "logl"),
trans=tmat2, M=5000)
wei1gom2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "gom", "logl"),
trans=tmat2, M=5000)
gom1gom2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "gom", "logl"),
trans=tmat2, M=5000)
gam1gom2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gom", "logl"),
trans=tmat2, M=5000)
logn1gom2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "gom", "logl"),
trans=tmat2, M=5000)
logl1gom2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "gom", "logl"),
trans=tmat2, M=5000)
##### gam2logl3
exp1gam2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "gam", "logl"),
trans=tmat2, M=5000)
wei1gam2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "gam", "logl"),
trans=tmat2, M=5000)
gom1gam2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "gam", "logl"),
trans=tmat2, M=5000)
gam1gam2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gam", "logl"),
trans=tmat2, M=5000)
logn1gam2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "gam", "logl"),
trans=tmat2, M=5000)
logl1gam2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "gam", "logl"),
trans=tmat2, M=5000)
##### logn2logl3
exp1logn2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "logn", "logl"),
trans=tmat2, M=5000)
wei1logn2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "logn", "logl"),
trans=tmat2, M=5000)
gom1logn2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "logn", "logl"),
trans=tmat2, M=5000)
gam1logn2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logn", "logl"),
trans=tmat2, M=5000)
logn1logn2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "logn", "logl"),
trans=tmat2, M=5000)
logl1logn2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "logn", "logl"),
trans=tmat2, M=5000)
##### logl2logl3
exp1logl2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "logl", "logl"),
trans=tmat2, M=5000)
wei1logl2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "logl", "logl"),
trans=tmat2, M=5000)
gom1logl2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "logl", "logl"),
trans=tmat2, M=5000)
gam1logl2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logl", "logl"),
trans=tmat2, M=5000)
logn1logl2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "logl", "logl"),
trans=tmat2, M=5000)
logl1logl2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "logl", "logl"),
trans=tmat2, M=5000)
##############################
### logn3
#############################
##### exp2logn3
exp1exp2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "exp", "logn"),
trans=tmat2, M=5000)
wei1exp2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "exp", "logn"),
trans=tmat2, M=5000)
gom1exp2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "exp", "logn"),
trans=tmat2, M=5000)
gam1exp2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "exp", "logn"),
trans=tmat2, M=5000)
logn1exp2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "exp", "logn"),
trans=tmat2, M=5000)
logl1exp2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "exp", "logn"),
trans=tmat2, M=5000)
##### wei2logn3
exp1wei2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "wei", "logn"),
trans=tmat2, M=5000)
wei1wei2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "wei", "logn"),
trans=tmat2, M=5000)
gom1wei2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "wei", "logn"),
trans=tmat2, M=5000)
gam1wei2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "wei", "logn"),
trans=tmat2, M=5000)
logn1wei2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "wei", "logn"),
trans=tmat2, M=5000)
logl1wei2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "wei", "logn"),
trans=tmat2, M=5000)
##### gom2logn3
exp1gom2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "gom", "logn"),
trans=tmat2, M=5000)
wei1gom2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "gom", "logn"),
trans=tmat2, M=5000)
gom1gom2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "gom", "logn"),
trans=tmat2, M=5000)
gam1gom2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gom", "logn"),
trans=tmat2, M=5000)
logn1gom2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "gom", "logn"),
trans=tmat2, M=5000)
logl1gom2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "gom", "logn"),
trans=tmat2, M=5000)
##### gam2logn3
exp1gam2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "gam", "logn"),
trans=tmat2, M=5000)
wei1gam2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "gam", "logn"),
trans=tmat2, M=5000)
gom1gam2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "gam", "logn"),
trans=tmat2, M=5000)
gam1gam2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gam", "logn"),
trans=tmat2, M=5000)
logn1gam2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "gam", "logn"),
trans=tmat2, M=5000)
logl1gam2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "gam", "logn"),
trans=tmat2, M=5000)
##### logn2logn3
exp1logn2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "logn", "logn"),
trans=tmat2, M=5000)
wei1logn2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "logn", "logn"),
trans=tmat2, M=5000)
gom1logn2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "logn", "logn"),
trans=tmat2, M=5000)
gam1logn2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logn", "logn"),
trans=tmat2, M=5000)
logn1logn2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "logn", "logn"),
trans=tmat2, M=5000)
logl1logn2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "logn", "logn"),
trans=tmat2, M=5000)
##### logl2logn3
exp1logl2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "logl", "logn"),
trans=tmat2, M=5000)
wei1logl2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "logl", "logn"),
trans=tmat2, M=5000)
gom1logl2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "logl", "logn"),
trans=tmat2, M=5000)
gam1logl2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logl", "logn"),
trans=tmat2, M=5000)
logn1logl2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "logl", "logn"),
trans=tmat2, M=5000)
logl1logl2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "logl", "logn"),
trans=tmat2, M=5000)
HealthEconomicsHackathon/ceatutorial documentation built on Nov. 7, 2019, 12:09 a.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
##########################################################################################
##========================================================================================
## CALCULATING STATE OCCUPANCY PROBABILITIES
##========================================================================================
##########################################################################################
### DUE TO COMPUTATIONAL DIFFICULTIES USING GOMPERTZ AND GENERALISED GAMMA FOR THE FIRST
### TRANSITION, MODELS THAT INCLUDE THIS USE MORE TIME POINTS. ALL OTHER MODELS USE THE
### EQUIVALENT OF MONTHLY CYCLES.
##############################
### RFC TREATMENT ARM
#############################
##############################
### exp3
#############################
ptm <- proc.time() #start counting the time it takes
##### exp2exp3
exp1exp2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "exp", "exp"),
trans=tmat2, M=5000)
(proc.time() - ptm) #stop counting the time it took (seconds)
(proc.time() - ptm)/60 #stop counting the time it took (minutes)
(proc.time() - ptm)/3600 #stop counting the time it took (hours)
wei1exp2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "exp", "exp"),
trans=tmat2, M=5000)
gom1exp2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
dist=cbind("gom", "exp", "exp"),
trans=tmat2, M=5000)
gam1exp2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "exp", "exp"),
trans=tmat2, M=5000)
logn1exp2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "exp", "exp"),
trans=tmat2, M=5000)
logl1exp2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "exp", "exp"),
trans=tmat2, M=5000)
##### wei2exp3
exp1wei2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "wei", "exp"),
trans=tmat2, M=5000)
wei1wei2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "wei", "exp"),
trans=tmat2, M=5000)
gom1wei2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
dist=cbind("gom", "wei", "exp"),
trans=tmat2, M=5000)
gam1wei2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "wei", "exp"),
trans=tmat2, M=5000)
logn1wei2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "wei", "exp"),
trans=tmat2, M=5000)
logl1wei2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "wei", "exp"),
trans=tmat2, M=5000)
##### gom2exp3
exp1gom2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "gom", "exp"),
trans=tmat2, M=5000)
wei1gom2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "gom", "exp"),
trans=tmat2, M=5000)
gom1gom2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
dist=cbind("gom", "gom", "exp"),
trans=tmat2, M=5000)
gam1gom2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gom", "exp"),
trans=tmat2, M=5000)
logn1gom2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "gom", "exp"),
trans=tmat2, M=5000)
logl1gom2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "gom", "exp"),
trans=tmat2, M=5000)
##### gam2exp3
exp1gam2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "gam", "exp"),
trans=tmat2, M=5000)
wei1gam2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "gam", "exp"),
trans=tmat2, M=5000)
gom1gam2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
dist=cbind("gom", "gam", "exp"),
trans=tmat2, M=5000)
gam1gam2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gam", "exp"),
trans=tmat2, M=5000)
logn1gam2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "gam", "exp"),
trans=tmat2, M=5000)
logl1gam2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "gam", "exp"),
trans=tmat2, M=5000)
##### logn2exp3
exp1logn2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "logn", "exp"),
trans=tmat2, M=5000)
wei1logn2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "logn", "exp"),
trans=tmat2, M=5000)
gom1logn2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
dist=cbind("gom", "logn", "exp"),
trans=tmat2, M=5000)
gam1logn2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logn", "exp"),
trans=tmat2, M=5000)
logn1logn2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "logn", "exp"),
trans=tmat2, M=5000)
logl1logn2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "logn", "exp"),
trans=tmat2, M=5000)
##### logl2exp3
exp1logl2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "logl", "exp"),
trans=tmat2, M=5000)
wei1logl2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "logl", "exp"),
trans=tmat2, M=5000)
gom1logl2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
dist=cbind("gom", "logl", "exp"),
trans=tmat2, M=5000)
gam1logl2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logl", "exp"),
trans=tmat2, M=5000)
logn1logl2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "logl", "exp"),
trans=tmat2, M=5000)
logl1logl2exp3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "logl", "exp"),
trans=tmat2, M=5000)
##############################
### wei3
#############################
##### exp2wei3
exp1exp2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "exp", "wei"),
trans=tmat2, M=5000)
wei1exp2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "exp", "wei"),
trans=tmat2, M=5000)
gom1exp2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("gom", "exp", "wei"),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
trans=tmat2, M=5000)
gam1exp2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "exp", "wei"),
trans=tmat2, M=5000)
logn1exp2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "exp", "wei"),
trans=tmat2, M=5000)
logl1exp2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "exp", "wei"),
trans=tmat2, M=5000)
##### wei2wei3
exp1wei2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "wei", "wei"),
trans=tmat2, M=5000)
wei1wei2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "wei", "wei"),
trans=tmat2, M=5000)
gom1wei2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
dist=cbind("gom", "wei", "wei"),
trans=tmat2, M=5000)
gam1wei2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "wei", "wei"),
trans=tmat2, M=5000)
logn1wei2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "wei", "wei"),
trans=tmat2, M=5000)
logl1wei2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "wei", "wei"),
trans=tmat2, M=5000)
##### gom2wei3
exp1gom2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "gom", "wei"),
trans=tmat2, M=5000)
wei1gom2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "gom", "wei"),
trans=tmat2, M=5000)
gom1gom2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
dist=cbind("gom", "gom", "wei"),
trans=tmat2, M=5000)
gam1gom2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gom", "wei"),
trans=tmat2, M=5000)
logn1gom2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "gom", "wei"),
trans=tmat2, M=5000)
logl1gom2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "gom", "wei"),
trans=tmat2, M=5000)
##### gam2wei3
exp1gam2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "gam", "wei"),
trans=tmat2, M=5000)
wei1gam2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "gam", "wei"),
trans=tmat2, M=5000)
gom1gam2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
dist=cbind("gom", "gam", "wei"),
trans=tmat2, M=5000)
gam1gam2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gam", "wei"),
trans=tmat2, M=5000)
logn1gam2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "gam", "wei"),
trans=tmat2, M=5000)
logl1gam2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "gam", "wei"),
trans=tmat2, M=5000)
##### logn2wei3
exp1logn2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "logn", "wei"),
trans=tmat2, M=5000)
wei1logn2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "logn", "wei"),
trans=tmat2, M=5000)
gom1logn2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
dist=cbind("gom", "logn", "wei"),
trans=tmat2, M=5000)
gam1logn2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logn", "wei"),
trans=tmat2, M=5000)
logn1logn2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "logn", "wei"),
trans=tmat2, M=5000)
logl1logn2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "logn", "wei"),
trans=tmat2, M=5000)
##### logl2wei3
exp1logl2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "logl", "wei"),
trans=tmat2, M=5000)
wei1logl2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "logl", "wei"),
trans=tmat2, M=5000)
gom1logl2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
dist=cbind("gom", "logl", "wei"),
trans=tmat2, M=5000)
gam1logl2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logl", "wei"),
trans=tmat2, M=5000)
logn1logl2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "logl", "wei"),
trans=tmat2, M=5000)
logl1logl2wei3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "logl", "wei"),
trans=tmat2, M=5000)
##############################
### gom3
#############################
##### exp2gom3
exp1exp2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "exp", "gom"),
trans=tmat2, M=5000)
wei1exp2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "exp", "gom"),
trans=tmat2, M=5000)
gom1exp2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
dist=cbind("gom", "exp", "gom"),
trans=tmat2, M=5000)
gam1exp2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "exp", "gom"),
trans=tmat2, M=5000)
logn1exp2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "exp", "gom"),
trans=tmat2, M=5000)
logl1exp2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "exp", "gom"),
trans=tmat2, M=5000)
##### wei2gom3
exp1wei2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "wei", "gom"),
trans=tmat2, M=5000)
wei1wei2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "wei", "gom"),
trans=tmat2, M=5000)
gom1wei2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
dist=cbind("gom", "wei", "gom"),
trans=tmat2, M=5000)
gam1wei2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "wei", "gom"),
trans=tmat2, M=5000)
logn1wei2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "wei", "gom"),
trans=tmat2, M=5000)
logl1wei2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "wei", "gom"),
trans=tmat2, M=5000)
##### gom2gom3
exp1gom2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "gom", "gom"),
trans=tmat2, M=5000)
wei1gom2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "gom", "gom"),
trans=tmat2, M=5000)
gom1gom2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
dist=cbind("gom", "gom", "gom"),
trans=tmat2, M=5000)
gam1gom2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gom", "gom"),
trans=tmat2, M=5000)
logn1gom2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "gom", "gom"),
trans=tmat2, M=5000)
logl1gom2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "gom", "gom"),
trans=tmat2, M=5000)
##### gam2gom3
exp1gam2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "gam", "gom"),
trans=tmat2, M=5000)
wei1gam2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "gam", "gom"),
trans=tmat2, M=5000)
gom1gam2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
dist=cbind("gom", "gam", "gom"),
trans=tmat2, M=5000)
gam1gam2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gam", "gom"),
trans=tmat2, M=5000)
logn1gam2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "gam", "gom"),
trans=tmat2, M=5000)
logl1gam2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "gam", "gom"),
trans=tmat2, M=5000)
##### logn2gom3
exp1logn2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "logn", "gom"),
trans=tmat2, M=5000)
wei1logn2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "logn", "gom"),
trans=tmat2, M=5000)
gom1logn2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
dist=cbind("gom", "logn", "gom"),
trans=tmat2, M=5000)
gam1logn2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logn", "gom"),
trans=tmat2, M=5000)
logn1logn2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "logn", "gom"),
trans=tmat2, M=5000)
logl1logn2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "logn", "gom"),
trans=tmat2, M=5000)
##### logl2gom3
exp1logl2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "logl", "gom"),
trans=tmat2, M=5000)
wei1logl2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "logl", "gom"),
trans=tmat2, M=5000)
gom1logl2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
dist=cbind("gom", "logl", "gom"),
trans=tmat2, M=5000)
gam1logl2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logl", "gom"),
trans=tmat2, M=5000)
logn1logl2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "logl", "gom"),
trans=tmat2, M=5000)
logl1logl2gom3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "logl", "gom"),
trans=tmat2, M=5000)
##############################
### gam3
#############################
##### exp2gam3
exp1exp2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "exp", "gam"),
trans=tmat2, M=5000)
wei1exp2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "exp", "gam"),
trans=tmat2, M=5000)
gom1exp2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("gom", "exp", "gam"),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
trans=tmat2, M=5000)
gam1exp2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "exp", "gam"),
trans=tmat2, M=5000)
logn1exp2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "exp", "gam"),
trans=tmat2, M=5000)
logl1exp2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "exp", "gam"),
trans=tmat2, M=5000)
##### wei2gam3
exp1wei2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "wei", "gam"),
trans=tmat2, M=5000)
wei1wei2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "wei", "gam"),
trans=tmat2, M=5000)
gom1wei2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("gom", "wei", "gam"),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
trans=tmat2, M=5000)
gam1wei2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "wei", "gam"),
trans=tmat2, M=5000)
logn1wei2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "wei", "gam"),
trans=tmat2, M=5000)
logl1wei2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "wei", "gam"),
trans=tmat2, M=5000)
##### gom2gam3
exp1gom2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "gom", "gam"),
trans=tmat2, M=5000)
wei1gom2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "gom", "gam"),
trans=tmat2, M=5000)
gom1gom2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("gom", "gom", "gam"),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
trans=tmat2, M=5000)
gam1gom2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gom", "gam"),
trans=tmat2, M=5000)
logn1gom2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "gom", "gam"),
trans=tmat2, M=5000)
logl1gom2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "gom", "gam"),
trans=tmat2, M=5000)
##### gam2gam3
exp1gam2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "gam", "gam"),
trans=tmat2, M=5000)
wei1gam2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "gam", "gam"),
trans=tmat2, M=5000)
gom1gam2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("gom", "gam", "gam"),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
trans=tmat2, M=5000)
gam1gam2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gam", "gam"),
trans=tmat2, M=5000)
logn1gam2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "gam", "gam"),
trans=tmat2, M=5000)
logl1gam2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "gam", "gam"),
trans=tmat2, M=5000)
##### logn2gam3
exp1logn2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "logn", "gam"),
trans=tmat2, M=5000)
wei1logn2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "logn", "gam"),
trans=tmat2, M=5000)
gom1logn2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("gom", "logn", "gam"),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
trans=tmat2, M=5000)
gam1logn2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logn", "gam"),
trans=tmat2, M=5000)
logn1logn2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "logn", "gam"),
trans=tmat2, M=5000)
logl1logn2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "logn", "gam"),
trans=tmat2, M=5000)
##### logl2gam3
exp1logl2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "logl", "gam"),
trans=tmat2, M=5000)
wei1logl2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "logl", "gam"),
trans=tmat2, M=5000)
gom1logl2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("gom", "logl", "gam"),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
trans=tmat2, M=5000)
gam1logl2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logl", "gam"),
trans=tmat2, M=5000)
logn1logl2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "logl", "gam"),
trans=tmat2, M=5000)
logl1logl2gam3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "logl", "gam"),
trans=tmat2, M=5000)
##############################
### logl3
#############################
##### exp2logl3
exp1exp2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "exp", "logl"),
trans=tmat2, M=5000)
wei1exp2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "exp", "logl"),
trans=tmat2, M=5000)
gom1exp2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("gom", "exp", "logl"),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
trans=tmat2, M=5000)
gam1exp2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "exp", "logl"),
trans=tmat2, M=5000)
logn1exp2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "exp", "logl"),
trans=tmat2, M=5000)
logl1exp2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "exp", "logl"),
trans=tmat2, M=5000)
##### wei2logl3
exp1wei2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "wei", "logl"),
trans=tmat2, M=5000)
wei1wei2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "wei", "logl"),
trans=tmat2, M=5000)
gom1wei2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("gom", "wei", "logl"),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
trans=tmat2, M=5000)
gam1wei2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "wei", "logl"),
trans=tmat2, M=5000)
logn1wei2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "wei", "logl"),
trans=tmat2, M=5000)
logl1wei2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "wei", "logl"),
trans=tmat2, M=5000)
##### gom2logl3
exp1gom2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "gom", "logl"),
trans=tmat2, M=5000)
wei1gom2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "gom", "logl"),
trans=tmat2, M=5000)
gom1gom2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("gom", "gom", "logl"),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
trans=tmat2, M=5000)
gam1gom2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gom", "logl"),
trans=tmat2, M=5000)
logn1gom2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "gom", "logl"),
trans=tmat2, M=5000)
logl1gom2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "gom", "logl"),
trans=tmat2, M=5000)
##### gam2logl3
exp1gam2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "gam", "logl"),
trans=tmat2, M=5000)
wei1gam2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "gam", "logl"),
trans=tmat2, M=5000)
gom1gam2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("gom", "gam", "logl"),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
trans=tmat2, M=5000)
gam1gam2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gam", "logl"),
trans=tmat2, M=5000)
logn1gam2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "gam", "logl"),
trans=tmat2, M=5000)
logl1gam2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "gam", "logl"),
trans=tmat2, M=5000)
##### logn2logl3
exp1logn2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "logn", "logl"),
trans=tmat2, M=5000)
wei1logn2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "logn", "logl"),
trans=tmat2, M=5000)
gom1logn2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("gom", "logn", "logl"),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
trans=tmat2, M=5000)
gam1logn2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logn", "logl"),
trans=tmat2, M=5000)
logn1logn2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "logn", "logl"),
trans=tmat2, M=5000)
logl1logn2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "logn", "logl"),
trans=tmat2, M=5000)
##### logl2logl3
exp1logl2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "logl", "logl"),
trans=tmat2, M=5000)
wei1logl2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "logl", "logl"),
trans=tmat2, M=5000)
gom1logl2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("gom", "logl", "logl"),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
trans=tmat2, M=5000)
gam1logl2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logl", "logl"),
trans=tmat2, M=5000)
logn1logl2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "logl", "logl"),
trans=tmat2, M=5000)
logl1logl2logl3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "logl", "logl"),
trans=tmat2, M=5000)
##############################
### logn3
#############################
##### exp2logn3
exp1exp2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "exp", "logn"),
trans=tmat2, M=5000)
wei1exp2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "exp", "logn"),
trans=tmat2, M=5000)
gom1exp2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("gom", "exp", "logn"),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
trans=tmat2, M=5000)
gam1exp2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "exp", "logn"),
trans=tmat2, M=5000)
logn1exp2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "exp", "logn"),
trans=tmat2, M=5000)
logl1exp2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "exp", "logn"),
trans=tmat2, M=5000)
##### wei2logn3
exp1wei2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "wei", "logn"),
trans=tmat2, M=5000)
wei1wei2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "wei", "logn"),
trans=tmat2, M=5000)
gom1wei2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("gom", "wei", "logn"),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
trans=tmat2, M=5000)
gam1wei2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "wei", "logn"),
trans=tmat2, M=5000)
logn1wei2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "wei", "logn"),
trans=tmat2, M=5000)
logl1wei2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "wei", "logn"),
trans=tmat2, M=5000)
##### gom2logn3
exp1gom2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "gom", "logn"),
trans=tmat2, M=5000)
wei1gom2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "gom", "logn"),
trans=tmat2, M=5000)
gom1gom2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("gom", "gom", "logn"),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
trans=tmat2, M=5000)
gam1gom2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gom", "logn"),
trans=tmat2, M=5000)
logn1gom2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "gom", "logn"),
trans=tmat2, M=5000)
logl1gom2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "gom", "logn"),
trans=tmat2, M=5000)
##### gam2logn3
exp1gam2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "gam", "logn"),
trans=tmat2, M=5000)
wei1gam2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "gam", "logn"),
trans=tmat2, M=5000)
gom1gam2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("gom", "gam", "logn"),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
trans=tmat2, M=5000)
gam1gam2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gam", "logn"),
trans=tmat2, M=5000)
logn1gam2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "gam", "logn"),
trans=tmat2, M=5000)
logl1gam2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "gam", "logn"),
trans=tmat2, M=5000)
##### logn2logn3
exp1logn2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "logn", "logn"),
trans=tmat2, M=5000)
wei1logn2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "logn", "logn"),
trans=tmat2, M=5000)
gom1logn2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("gom", "logn", "logn"),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
trans=tmat2, M=5000)
gam1logn2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logn", "logn"),
trans=tmat2, M=5000)
logn1logn2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "logn", "logn"),
trans=tmat2, M=5000)
logl1logn2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "logn", "logn"),
trans=tmat2, M=5000)
##### logl2logn3
exp1logl2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("exp", "logl", "logn"),
trans=tmat2, M=5000)
wei1logl2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("wei", "logl", "logn"),
trans=tmat2, M=5000)
gom1logl2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("gom", "logl", "logn"),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)),
trans=tmat2, M=5000)
gam1logl2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logl", "logn"),
trans=tmat2, M=5000)
logn1logl2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logn", "logl", "logn"),
trans=tmat2, M=5000)
logl1logl2logn3simRFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(1,1,1),
dist=cbind("logl", "logl", "logn"),
trans=tmat2, M=5000)
##############################
# ### FC TREATMENT ARM
# #############################
#
# ##############################
# ### exp3
# #############################
#
##### exp2exp3
exp1exp2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "exp", "exp"),
trans=tmat2, M=5000)
wei1exp2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "exp", "exp"),
trans=tmat2, M=5000)
gom1exp2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "exp", "exp"),
trans=tmat2, M=5000)
gam1exp2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)),
dist=cbind("gam", "exp", "exp"),
trans=tmat2, M=5000)
logn1exp2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "exp", "exp"),
trans=tmat2, M=5000)
logl1exp2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "exp", "exp"),
trans=tmat2, M=5000)
##### wei2exp3
exp1wei2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "wei", "exp"),
trans=tmat2, M=5000)
wei1wei2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "wei", "exp"),
trans=tmat2, M=5000)
gom1wei2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "wei", "exp"),
trans=tmat2, M=5000)
gam1wei2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "wei", "exp"),
trans=tmat2, M=5000)
logn1wei2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "wei", "exp"),
trans=tmat2, M=5000)
logl1wei2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "wei", "exp"),
trans=tmat2, M=5000)
##### gom2exp3
exp1gom2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "gom", "exp"),
trans=tmat2, M=5000)
wei1gom2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "gom", "exp"),
trans=tmat2, M=5000)
gom1gom2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "gom", "exp"),
trans=tmat2, M=5000)
gam1gom2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gom", "exp"),
trans=tmat2, M=5000)
logn1gom2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "gom", "exp"),
trans=tmat2, M=5000)
logl1gom2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "gom", "exp"),
trans=tmat2, M=5000)
##### gam2exp3
exp1gam2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "gam", "exp"),
trans=tmat2, M=5000)
wei1gam2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "gam", "exp"),
trans=tmat2, M=5000)
gom1gam2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "gam", "exp"),
trans=tmat2, M=5000)
gam1gam2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gam", "exp"),
trans=tmat2, M=5000)
logn1gam2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "gam", "exp"),
trans=tmat2, M=5000)
logl1gam2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "gam", "exp"),
trans=tmat2, M=5000)
##### logn2exp3
exp1logn2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "logn", "exp"),
trans=tmat2, M=5000)
wei1logn2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "logn", "exp"),
trans=tmat2, M=5000)
gom1logn2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "logn", "exp"),
trans=tmat2, M=5000)
gam1logn2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logn", "exp"),
trans=tmat2, M=5000)
logn1logn2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "logn", "exp"),
trans=tmat2, M=5000)
logl1logn2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "logn", "exp"),
trans=tmat2, M=5000)
##### logl2exp3
exp1logl2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "logl", "exp"),
trans=tmat2, M=5000)
wei1logl2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "logl", "exp"),
trans=tmat2, M=5000)
gom1logl2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "logl", "exp"),
trans=tmat2, M=5000)
gam1logl2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logl", "exp"),
trans=tmat2, M=5000)
logn1logl2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "logl", "exp"),
trans=tmat2, M=5000)
logl1logl2exp3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "logl", "exp"),
trans=tmat2, M=5000)
##############################
### wei3
#############################
##### exp2wei3
exp1exp2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "exp", "wei"),
trans=tmat2, M=5000)
wei1exp2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "exp", "wei"),
trans=tmat2, M=5000)
gom1exp2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "exp", "wei"),
trans=tmat2, M=5000)
gam1exp2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "exp", "wei"),
trans=tmat2, M=5000)
logn1exp2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "exp", "wei"),
trans=tmat2, M=5000)
logl1exp2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "exp", "wei"),
trans=tmat2, M=5000)
##### wei2wei3
exp1wei2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "wei", "wei"),
trans=tmat2, M=5000)
wei1wei2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "wei", "wei"),
trans=tmat2, M=5000)
gom1wei2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "wei", "wei"),
trans=tmat2, M=5000)
gam1wei2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "wei", "wei"),
trans=tmat2, M=5000)
logn1wei2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "wei", "wei"),
trans=tmat2, M=5000)
logl1wei2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "wei", "wei"),
trans=tmat2, M=5000)
##### gom2wei3
exp1gom2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "gom", "wei"),
trans=tmat2, M=5000)
wei1gom2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "gom", "wei"),
trans=tmat2, M=5000)
gom1gom2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "gom", "wei"),
trans=tmat2, M=5000)
gam1gom2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gom", "wei"),
trans=tmat2, M=5000)
logn1gom2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "gom", "wei"),
trans=tmat2, M=5000)
logl1gom2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "gom", "wei"),
trans=tmat2, M=5000)
##### gam2wei3
exp1gam2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "gam", "wei"),
trans=tmat2, M=5000)
wei1gam2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "gam", "wei"),
trans=tmat2, M=5000)
gom1gam2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "gam", "wei"),
trans=tmat2, M=5000)
gam1gam2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gam", "wei"),
trans=tmat2, M=5000)
logn1gam2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "gam", "wei"),
trans=tmat2, M=5000)
logl1gam2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "gam", "wei"),
trans=tmat2, M=5000)
##### logn2wei3
exp1logn2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "logn", "wei"),
trans=tmat2, M=5000)
wei1logn2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "logn", "wei"),
trans=tmat2, M=5000)
gom1logn2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "logn", "wei"),
trans=tmat2, M=5000)
gam1logn2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logn", "wei"),
trans=tmat2, M=5000)
logn1logn2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "logn", "wei"),
trans=tmat2, M=5000)
logl1logn2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "logn", "wei"),
trans=tmat2, M=5000)
##### logl2wei3
exp1logl2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "logl", "wei"),
trans=tmat2, M=5000)
wei1logl2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "logl", "wei"),
trans=tmat2, M=5000)
gom1logl2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "logl", "wei"),
trans=tmat2, M=5000)
gam1logl2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logl", "wei"),
trans=tmat2, M=5000)
logn1logl2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "logl", "wei"),
trans=tmat2, M=5000)
logl1logl2wei3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "logl", "wei"),
trans=tmat2, M=5000)
##############################
### gom3
#############################
##### exp2gom3
exp1exp2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "exp", "gom"),
trans=tmat2, M=5000)
wei1exp2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "exp", "gom"),
trans=tmat2, M=5000)
gom1exp2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "exp", "gom"),
trans=tmat2, M=5000)
gam1exp2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "exp", "gom"),
trans=tmat2, M=5000)
logn1exp2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "exp", "gom"),
trans=tmat2, M=5000)
logl1exp2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "exp", "gom"),
trans=tmat2, M=5000)
##### wei2gom3
exp1wei2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "wei", "gom"),
trans=tmat2, M=5000)
wei1wei2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "wei", "gom"),
trans=tmat2, M=5000)
gom1wei2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "wei", "gom"),
trans=tmat2, M=5000)
gam1wei2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "wei", "gom"),
trans=tmat2, M=5000)
logn1wei2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "wei", "gom"),
trans=tmat2, M=5000)
logl1wei2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "wei", "gom"),
trans=tmat2, M=5000)
##### gom2gom3
exp1gom2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "gom", "gom"),
trans=tmat2, M=5000)
wei1gom2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "gom", "gom"),
trans=tmat2, M=5000)
gom1gom2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "gom", "gom"),
trans=tmat2, M=5000)
gam1gom2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gom", "gom"),
trans=tmat2, M=5000)
logn1gom2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "gom", "gom"),
trans=tmat2, M=5000)
logl1gom2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "gom", "gom"),
trans=tmat2, M=5000)
##### gam2gom3
exp1gam2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "gam", "gom"),
trans=tmat2, M=5000)
wei1gam2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "gam", "gom"),
trans=tmat2, M=5000)
gom1gam2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "gam", "gom"),
trans=tmat2, M=5000)
gam1gam2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gam", "gom"),
trans=tmat2, M=5000)
logn1gam2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "gam", "gom"),
trans=tmat2, M=5000)
logl1gam2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "gam", "gom"),
trans=tmat2, M=5000)
##### logn2gom3
exp1logn2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "logn", "gom"),
trans=tmat2, M=5000)
wei1logn2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "logn", "gom"),
trans=tmat2, M=5000)
gom1logn2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "logn", "gom"),
trans=tmat2, M=5000)
gam1logn2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logn", "gom"),
trans=tmat2, M=5000)
logn1logn2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "logn", "gom"),
trans=tmat2, M=5000)
logl1logn2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "logn", "gom"),
trans=tmat2, M=5000)
##### logl2gom3
exp1logl2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "logl", "gom"),
trans=tmat2, M=5000)
wei1logl2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "logl", "gom"),
trans=tmat2, M=5000)
gom1logl2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "logl", "gom"),
trans=tmat2, M=5000)
gam1logl2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logl", "gom"),
trans=tmat2, M=5000)
logn1logl2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "logl", "gom"),
trans=tmat2, M=5000)
logl1logl2gom3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "logl", "gom"),
trans=tmat2, M=5000)
##############################
### gam3
#############################
##### exp2gam3
exp1exp2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "exp", "gam"),
trans=tmat2, M=5000)
wei1exp2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "exp", "gam"),
trans=tmat2, M=5000)
gom1exp2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "exp", "gam"),
trans=tmat2, M=5000)
gam1exp2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "exp", "gam"),
trans=tmat2, M=5000)
logn1exp2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "exp", "gam"),
trans=tmat2, M=5000)
logl1exp2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "exp", "gam"),
trans=tmat2, M=5000)
##### wei2gam3
exp1wei2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "wei", "gam"),
trans=tmat2, M=5000)
wei1wei2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "wei", "gam"),
trans=tmat2, M=5000)
gom1wei2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "wei", "gam"),
trans=tmat2, M=5000)
gam1wei2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "wei", "gam"),
trans=tmat2, M=5000)
logn1wei2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "wei", "gam"),
trans=tmat2, M=5000)
logl1wei2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "wei", "gam"),
trans=tmat2, M=5000)
##### gom2gam3
exp1gom2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "gom", "gam"),
trans=tmat2, M=5000)
wei1gom2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "gom", "gam"),
trans=tmat2, M=5000)
gom1gom2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "gom", "gam"),
trans=tmat2, M=5000)
gam1gom2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gom", "gam"),
trans=tmat2, M=5000)
logn1gom2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "gom", "gam"),
trans=tmat2, M=5000)
logl1gom2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "gom", "gam"),
trans=tmat2, M=5000)
##### gam2gam3
exp1gam2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "gam", "gam"),
trans=tmat2, M=5000)
wei1gam2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "gam", "gam"),
trans=tmat2, M=5000)
gom1gam2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "gam", "gam"),
trans=tmat2, M=5000)
gam1gam2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gam", "gam"),
trans=tmat2, M=5000)
logn1gam2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "gam", "gam"),
trans=tmat2, M=5000)
logl1gam2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "gam", "gam"),
trans=tmat2, M=5000)
##### logn2gam3
exp1logn2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "logn", "gam"),
trans=tmat2, M=5000)
wei1logn2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "logn", "gam"),
trans=tmat2, M=5000)
gom1logn2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "logn", "gam"),
trans=tmat2, M=5000)
gam1logn2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logn", "gam"),
trans=tmat2, M=5000)
logn1logn2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "logn", "gam"),
trans=tmat2, M=5000)
logl1logn2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "logn", "gam"),
trans=tmat2, M=5000)
##### logl2gam3
exp1logl2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "logl", "gam"),
trans=tmat2, M=5000)
wei1logl2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "logl", "gam"),
trans=tmat2, M=5000)
gom1logl2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "logl", "gam"),
trans=tmat2, M=5000)
gam1logl2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logl", "gam"),
trans=tmat2, M=5000)
logn1logl2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "logl", "gam"),
trans=tmat2, M=5000)
logl1logl2gam3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "logl", "gam"),
trans=tmat2, M=5000)
##############################
### logl3
#############################
##### exp2logl3
exp1exp2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "exp", "logl"),
trans=tmat2, M=5000)
wei1exp2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "exp", "logl"),
trans=tmat2, M=5000)
gom1exp2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "exp", "logl"),
trans=tmat2, M=5000)
gam1exp2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "exp", "logl"),
trans=tmat2, M=5000)
logn1exp2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "exp", "logl"),
trans=tmat2, M=5000)
logl1exp2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "exp", "logl"),
trans=tmat2, M=5000)
##### wei2logl3
exp1wei2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "wei", "logl"),
trans=tmat2, M=5000)
wei1wei2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "wei", "logl"),
trans=tmat2, M=5000)
gom1wei2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "wei", "logl"),
trans=tmat2, M=5000)
gam1wei2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "wei", "logl"),
trans=tmat2, M=5000)
logn1wei2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "wei", "logl"),
trans=tmat2, M=5000)
logl1wei2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "wei", "logl"),
trans=tmat2, M=5000)
##### gom2logl3
exp1gom2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "gom", "logl"),
trans=tmat2, M=5000)
wei1gom2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "gom", "logl"),
trans=tmat2, M=5000)
gom1gom2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "gom", "logl"),
trans=tmat2, M=5000)
gam1gom2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gom", "logl"),
trans=tmat2, M=5000)
logn1gom2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "gom", "logl"),
trans=tmat2, M=5000)
logl1gom2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "gom", "logl"),
trans=tmat2, M=5000)
##### gam2logl3
exp1gam2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "gam", "logl"),
trans=tmat2, M=5000)
wei1gam2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "gam", "logl"),
trans=tmat2, M=5000)
gom1gam2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "gam", "logl"),
trans=tmat2, M=5000)
gam1gam2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gam", "logl"),
trans=tmat2, M=5000)
logn1gam2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "gam", "logl"),
trans=tmat2, M=5000)
logl1gam2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "gam", "logl"),
trans=tmat2, M=5000)
##### logn2logl3
exp1logn2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "logn", "logl"),
trans=tmat2, M=5000)
wei1logn2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "logn", "logl"),
trans=tmat2, M=5000)
gom1logn2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "logn", "logl"),
trans=tmat2, M=5000)
gam1logn2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logn", "logl"),
trans=tmat2, M=5000)
logn1logn2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "logn", "logl"),
trans=tmat2, M=5000)
logl1logn2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "logn", "logl"),
trans=tmat2, M=5000)
##### logl2logl3
exp1logl2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "logl", "logl"),
trans=tmat2, M=5000)
wei1logl2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "logl", "logl"),
trans=tmat2, M=5000)
gom1logl2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "logl", "logl"),
trans=tmat2, M=5000)
gam1logl2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logl", "logl"),
trans=tmat2, M=5000)
logn1logl2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "logl", "logl"),
trans=tmat2, M=5000)
logl1logl2logl3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "logl", "logl"),
trans=tmat2, M=5000)
##############################
### logn3
#############################
##### exp2logn3
exp1exp2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "exp", "logn"),
trans=tmat2, M=5000)
wei1exp2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "exp", "logn"),
trans=tmat2, M=5000)
gom1exp2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "exp", "logn"),
trans=tmat2, M=5000)
gam1exp2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "exp", "logn"),
trans=tmat2, M=5000)
logn1exp2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "exp", "logn"),
trans=tmat2, M=5000)
logl1exp2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "exp", "logn"),
trans=tmat2, M=5000)
##### wei2logn3
exp1wei2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "wei", "logn"),
trans=tmat2, M=5000)
wei1wei2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "wei", "logn"),
trans=tmat2, M=5000)
gom1wei2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "wei", "logn"),
trans=tmat2, M=5000)
gam1wei2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "wei", "logn"),
trans=tmat2, M=5000)
logn1wei2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "wei", "logn"),
trans=tmat2, M=5000)
logl1wei2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "wei", "logn"),
trans=tmat2, M=5000)
##### gom2logn3
exp1gom2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "gom", "logn"),
trans=tmat2, M=5000)
wei1gom2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "gom", "logn"),
trans=tmat2, M=5000)
gom1gom2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "gom", "logn"),
trans=tmat2, M=5000)
gam1gom2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gom", "logn"),
trans=tmat2, M=5000)
logn1gom2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "gom", "logn"),
trans=tmat2, M=5000)
logl1gom2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "gom", "logn"),
trans=tmat2, M=5000)
##### gam2logn3
exp1gam2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "gam", "logn"),
trans=tmat2, M=5000)
wei1gam2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "gam", "logn"),
trans=tmat2, M=5000)
gom1gam2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "gam", "logn"),
trans=tmat2, M=5000)
gam1gam2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "gam", "logn"),
trans=tmat2, M=5000)
logn1gam2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "gam", "logn"),
trans=tmat2, M=5000)
logl1gam2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "gam", "logn"),
trans=tmat2, M=5000)
##### logn2logn3
exp1logn2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "logn", "logn"),
trans=tmat2, M=5000)
wei1logn2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "logn", "logn"),
trans=tmat2, M=5000)
gom1logn2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "logn", "logn"),
trans=tmat2, M=5000)
gam1logn2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logn", "logn"),
trans=tmat2, M=5000)
logn1logn2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "logn", "logn"),
trans=tmat2, M=5000)
logl1logn2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "logn", "logn"),
trans=tmat2, M=5000)
##### logl2logn3
exp1logl2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("exp", "logl", "logn"),
trans=tmat2, M=5000)
wei1logl2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("wei", "logl", "logn"),
trans=tmat2, M=5000)
gom1logl2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,100/12,1/12), seq(100/12+1/144, 12, 1/144), seq(12+1/600,15, 1/600)), dist=cbind("gom", "logl", "logn"),
trans=tmat2, M=5000)
gam1logl2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
timeseq_ext=c(seq(49/12,144/12,1/12), seq(144/12+1/144, 15, 1/144)), dist=cbind("gam", "logl", "logn"),
trans=tmat2, M=5000)
logn1logl2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logn", "logl", "logn"),
trans=tmat2, M=5000)
logl1logl2logn3simFC<-semiMarkov(ncovs=c(1,1,1),covs=rbind("treat", "treat", "treat"),
coveval=rbind(0,0,0),
dist=cbind("logl", "logl", "logn"),
trans=tmat2, M=5000)
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.